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I ' I We analyse the constraints on dimension-six AF = 1 effective operators in models respect- 

^^ ■ ing the MFV hypothesis, both in the onc-Higgs doublet case and in the two-Higgs doublet 

00 I scenario with large tan/3. The constraints are derived mainly from the b ^ s inclusive ob- 

^D ■ servables measured at the B factories. The implications of these bounds in view of improved 

measurements in exclusive and inclusive observables in fe ^ s£^£~ and s -^ dvv transitions 
, , are discussed. 
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1 Introduction 

The Standard Model (SM) can be viewed as the renormahzable part of an effective field theory, 
valid up to some still undetermined cut-off scale A above the electroweak scale, {■\/2Gf)~^''^ ~ 
250 GeV. Theoretical arguments based on a natural solution of the hierarchy problem suggest 
that A should not exceed a few TeV. This expectation leads to a paradox when combined with 
the absence of significant deviations from the SM in loop-induced flavour-violating observables, 
potentially sensitive to very high energy scales. An effective solution to this problem is provided 
by the so-called hypothesis of Minimal Flavour Violation [1] , namely by the assumption that the 
SM Yukawa couplings are the only relevant breaking sources of the SU{'if' flavour symmetry [2] 
of the low-energy effective theory.^ 

This symmetry and symmetry-breaking ansatz is realised in various explicit extensions of the 
SM, such as supersymmetric models (see e.g. Ref [5,6] and [7]) or models with extra dimensions 
(see e.g. Ref. [8]). However, its main virtue is the possibility to perform a general analysis of 
new-physics effects in low-energy observables independently of the ultraviolet completion of the 
model. As shown in Ref. [1], the MFV hypothesis allows to build a rather predictive effective 
theory in terms of SM and Higgs fields. The predictions on flavour-violating observables derived 
within this effective theory are powerful tests of the underlying flavour structure of the model: if 
falsified, these tests would unambiguously signal the presence of new symmetry-breaking terms. 

The observables most relevant to test the MFV hypothesis and, within this framework, to 
constrain the structure of the effective theory are AF = 2 and A.F = 1 flavour-changing neutral- 
current (FCNC) processes. An updated analysis of the /S.F = 2 sector, or the meson-antimeson 
mixing amplitudes, has been presented recently in Ref. [9]. The goal of this work is a complete 
analysis of the AF = 1 sector, or the rare-decay amplitudes. 

Using the currently available measurements of AF = 1 FCNC processes from 6 — > s and 
s ^ d transitions (see Table 1) we derive updated bounds on the effective scale of new physics 
within MFV models. We consider in particular both the scenario of one effective Higgs doublet 
and the case of two Higgs doublets and large tan/3, where we are free to change the relative 
normalization of the two Yukawa couplings and to decouple the breaking of C/(1)pq and 5C/(3)^ 
global symmetries [1]. 

Having derived the bounds on the effective operators from the observables listed in Table 1, 
we derive a series of predictions for exclusive and inclusive observables in 6 — > s^~^l~ and s — > dvU 
transitions which have not been measured so far with high accuracy. On the one hand, these 
predictions indicate where to look for large new physics effects in the flavour sector, even under 
the pessimistic hypothesis of MFV. On the other hand, some of these predictions could provide, 
in the future, a proof of the MFV hypothesis: a set of deviations from the SM exhibiting the 
correlation predicted by this symmetry structure. 

The paper is organised as follows: in Section 2 we review the structure of the effective AF = 1 
Hamiltonian under the MFV hypothesis. The theoretical expressions of the observables analysed, 
in terms of the Wilson coefficients of the effective theory, are presented in Section 3. The numerical 
bounds on the scale of new physics and the predictions for future measurements are discussed in 
Section 4 and 5, respectively. The results are briefly summarised in the Conclusions. 



^ For earlier/alternative definitions of the MFV hypothesis see Ref. [3,4]. 



Observable Experiment SM prediction 

BjB ^ Xs7)iE-,>i.6 Gcvi (3.52 ± 0.24) x 10"^ [10] (3.13 ± 0.23) x 10"^ [11-13] 

S(i?^X,£+r)[^.g[o.04,i.o]GcV^] (0.6±0.5)xl0-« (0.8 ± 0.2) x 10"^ 

i3(5 ^ ^.^+^~)[g2g[i.o,6.o] Gcv2] (1-6 ± 0.5) X 10-6 [14, 15]'' (1.6 ± 0.1) x 10^^ [16-23] 

B{B -^ X,£+r )^g2>i4,4 Gev^] (4.4 ± 1.3) X 10"^ (2.4 ± 0.8) x 10"^ 

BjBs -^ fi+fj-) < 5.8 X 10"^ (95% CL) [24] (4.1 ± 0.8) x lO"'^ [25,26] 

ApBiB ^ ii'*£+r)[^2<6.25 Gev^] 0.24+^;^^ -0.01 ±0.02 

I^,,(i? ^ i^*^+^-)[g2>io.24 Gev^] 0-76;°:|l [27] 0.20 ± 0.08 [28-31] 

B{K+ -^ -K+uv) (14.7+8;!)^) X 10^11 [32] (8.6 ± 0.9) x 10"" [33-37] 



"Here we quote naive averages of the values obtained by the experiments and with symmetrized errors. 

Table 1: Main observables used to determine bounds on the MFV dimension-six operators. The 
SM predictions are updated according to the most recent determinations of the SM input values 
(see Section 4). 

2 MFV and AF = 1 processes 

Under the MFV hypothesis, the dimension-six effective operators relevant to down-type FCNC 
transitions, both with one or two Higgs doublets, can be defined as follows [1]: 

Ohi = i (QlAfcTmQl) h\jD^Hu , Oh2 = i (OlXfct'^j^Ql) H^t'^D^Hu , 

Ogi = Hd {QLXFcXda^uT-DR} (^.G^J , Og2 = {QlXyc1^T-Ql) {gsD^Gl,) , 

Ofi = Hd {QL\Fc\dCJ^uDR) (eF^^) , Of2 = {Ql>^fc1^Ql) {eD^F^^) , (1) 

Oil = {QlXfcI^Ql) (LlI^Ll) , Oe2 = (OlAfc7m^"Ql) (L^Tm^^^l) > 

Oi3 = {QlXfcI^Ql) (ER-f^ER) , Osi = {QL^FC^dDR) {ErX^Ll) . 

The five fermion fields {Ql, Dr, Ur, Ll, Er) carry a flavour index (i = 1, ... ,3): since we are 
interested in FCNC processes of down-type quarks, it is convenient to work in the mass-eigenstate 
basis of charged leptons and down-type quarks, where X^^e = diag{md//{H£))}- In this basis the 
flavour-changing coupling (AFc)ij can be expanded as 

(Apc) = |M)«-''''=''^''- '^'- (2) 

I i = j , 

in terms of the top-quark Yukawa coupling (Af = mt/{Hu) ~ 1) and the CKM matrix (V^j). Note 
that we have defined the operators linear in the gauge fields including appropriate powers of the 
corresponding gauge couplings (contrary to the original definition of Ref. [1]). Moreover, we have 
included from the beginning the scalar-density operator ©51, which plays a relevant role in the 
two-Higgs doublet case at large tan/3 = {Hjj) /{Hd)? For all the operators linear in A^ we could 
also consider a corresponding set with opposite chirality {Dr <-^ Ql)] however, the hierarchical 
structure of A^ implies that only one chirality structure is relevant. In the following we restrict 
the attention to flavour transitions of the type j > i, such as the leading chirality structure is the 
one in Eq. (1). 



^ In principle, in the two-Higgs doublet case we can consider additional operators obtained from those in (1) by 
the exchange Hu -^ Hd and/or Afc — * YdY^Xfc- However, for all the B physics observables we analyse in this 
work, the effects of these additional operators can be reabsorbed into a redefinition of the couplings of the operators 
in (1). 



Before the breaking of the electroweak symmetry, we define the MFV effective Lagrangian 
as jI^^Jmfy ~ iJ2n'^riOn)/^'^, where Oj are 0(1) couphngs and A is the effective scale of new 
physics. After the spontaneous breaking of the electroweak symmetry, the operators in (1) are 
mapped onto the standard basis of FCNC operators, defined by 



GpOc 

Itt sin" Ww 



Ksij>i] - - ^^^ ^;„2 fl... ^3i^3j 2^ C„ Qn + h.C. , (3) 



where 



Q? = jimj dia^iy (1 + 75) dj F^^ , Qg = prrij dia^^T"" (1 + 75) dj G"^^ , 

Qg = da^ (1 - 75) dj J2e hfi^ , Qio = da^ (1 - 75) dj Y,i h^ild , (4) 

Qpu = dilfi (1 - 75) dj J2u ^7m (1 - 75) i^ , Qs = di (1 + 75) dj e{l - 75) ^ • 

Defining 

-— ] Qi , Ao = ^ 2.4 TeV , (5) 

A J Oem 

the modified initial conditions of the Wilson coefficients of this effective Hamiltonian at the elec- 
troweak scale, 6Ci = Cf^"^ {im) - Cf^{m), are [1] 

5Cj = 2g\Fi , 5Cg = 2(7'eGi , 

dCg = en + ei>2 + e^s - [(1 - 4sin2 Ow^z + '2e'^eF2] , SCw = ez - en - ei2 + e^s , (6) 



■^j 



where ez = (e^i + eH2)/2. 

At this point it is useful to compare our effective AF = 1 Hamiltonian to those adopted in 
similar analyses in the recent literature. Bounds on the scale of new physics and corresponding 
predictions for rare decays in MFV scenarios have been discussed also by other authors (see 
e.g. Ref. [38-40]). However, most of the recent analyses concentrated on a specific version of the 
general MFV framework [1], the so-called constrained MFY (CMFV) [3]. While MFV as proposed 
in [1] is an hypothesis about the symmetry-breaking structure of the SU{3)^ flavour symmetry, the 
CMFV contains a further dynamical assumption: the hypothesis of no new effective dimension- 
six flavour-changing operators beside the SM ones (after electroweak symmetry breaking). In 
practice, also this dynamical assumption can be related to a symmetry-breaking issue: namely the 
hypothesis about the breaking of the [/(1)pq symmetry of the SM gauge Lagrangian [1]. Indeed 
independently from the value of tan/3 = {Hu)/{H[)), in absence of a sizable C/(1)pq breaking 
the coefficient of the scalar operator is too small to compensate the strong A^.A^ suppression in 
Eq. (6). Since large C/(1)pq breakings can occur in well-motivated extensions of the SM, such as 
supersymmetric models (see e.g. [41-43]), we make no specific assumptions about its size and take 
into account also the effects of the scalar operators. For the same reason, using the definition of 
Apc in (2) we remove from the operator basis terms contributing to flavour-conserving processes 
induced by the diagonal component of YijY^. At large tan/? flavour-diagonal terms could receive 
additional contributions of the type Yd^^h therefore we cannot use them to bound the flavour- 
changing terms. This does not occur in the CMFV, where Z ^ bb provide a useful constraint on 
flavour-changing operators [40]. 

2.1 Effective weak Hamiltonian at fi = rrih 

When computing AF = 1 observables, we need to take into account also the contributions of four- 
quark operators. Following the notation of Ref. [16], we write the complete effective Hamiltonian 



relevant to & — > s transitions as 
AC 



-^ \ Ec^^-^^f- + Ks^cfc)^f + y:sythC'^p^ + yLytKp'^ \ + h.c. , (?) 



where 



P^ = {slJ^T-ul){ulj''T%l) , P5 = (sl7/.i7/.2 7/.3^l) E,(97'^^7'^^7^^g) , 

P| = {sLl^UL){uLl^bL) , Pe = (sL7w7/.2 7/.3r"^i) E5(^7'^^7^^7^^T'''7) , 

Pf = {sl7^.T''cl){cl7''T%l) , Pi = jimb{sLa^'''hR)F^., , 

Pi = {sLl^.CL){cLl^bL) , Ps = j-mb{sL(T^''T''bR)G^,, , (8) 

Pa = {sLl>.T-bL) E,{Q7^T-q) , ^9 = ^(sL^/^bL) Eiih^^^) ■ 

With respect to the SM hterature we also add the scalar-density operator with right-handed b 
quark, 

P^ = T^{sLbR){eR£L) , (9) 

which plays an important role in the large tan /3 regime. This operator is present also in the SM, 
but it is usually neglected because of the of strong suppression of its Wilson coefficients. 

In principle we should consider new-physics effects both in the four-quark (Pi^q) and in the 
FCNC operators. However, four-quark operators receive a large SM contribution which is natu- 
rally much larger with respect to the new-physics one in the MFV scenario. Moreover four-quark 
operators do not contribute at the tree-level to the observables we are considering (FCNC pro- 
cesses). As a result, it is a good approximation to consider new-physics effects only in the leading 
FCNC operators:^ 

Ckf^H) = ClsM{m) + -r^^SCg, (10) 

47r sm t^w 

4tt sm t^w 

with the 6Ci given in Eq. (6). The initial conditions for the SM terms are obtained by a NNLO 
matching of the effective Wilson coefficients at the high scale fin G [80, 160] GeV, using CKM 
unitarity to get rid of V*i^Vcs/V^lVts [16]. 



^ Even including explicit new-physics contributions to Ci-e, to an excellent accuracy these could reabsorbed 
into the Wilson coefficients of the FCNC operators. 



Quantity Value 



sin^ 9w 


0.2312 


ocs{mz) 


0.1187 


aem{mz) 


1/127.8 



m^^inib) 4.2 GeV 
mj^{mt) 165 GeV 

Table 2: Reference input values used in the matching and in the RGE evolution of the Wilson 
coefficients. 



To a good approximation, in all the considered low-energy observables dCj and 5Cs appear in 
a fixed linear combination: 

6C7 + 0.3 SCs . (11) 

For this reason, in the following we set 5Cs = and treat only 50^ as independent fit parameter 
(avoiding the flat direction in the parameter space determined by Eq. (11)). The bounds on 6C7 
thus obtained should therefore be interpreted as bounds on the 6Cj-6Cs combination in Eq. (11). 
The Wilson coefficients are evolved from the high scale down to fib = 0{mf,) using SM renor- 
malization group equations at the NNLO accuracy [44,45]. At the low scale it is convenient to 
define the effective coefficients [18] 

A7 = -^,Cr{fib)-\c3{fib)-^C,{f,b)-^C,{f,b)-^Ce{fib), (12) 

as{fJ,bj 3 9 3 9 

A9 = ^^C9(/i,) + X]f^.(^f<)7j?ln^ + ^C73(/.fe) + ^C5(/i,) + ^C6(Mfe), (13) 
«s(/^6) ~^ A* 3 9 27 

47r 
^10 = -^^Cioifib), (14) 

Tg = ^Ci(/Zb)+C2(/i6) + 6C3(Mfe)+60C5(/ife), (15) 

7 2 32 

U9 = -^C3{f,b)-^C^if^b)-38C5ifib)-Y'^6if^b), (16) 

1 2 32 

Wg = --C3(A^6)-3C4(A^6)-8C5(Mfe)-yC6(/Xb), (17) 

where Ci{fi) = Cj{iJ,) — Cf{iJ,), j^g is given in [16] and we have neglected the tiny u-quark loop 
contribution to Ag. The reference input values used in this procedure are collected in Table 2. 

Beside B physics, we are interested also in s — > d transitions and particularly in the rare 
decays K -^ ttvu. In the MFV framework these are described by the following simple effective 
Hamiltonian 

t=e,fj.,T 

where 

Puu = is-f^d){Pa^il--f5)i^i) , (19) 

y. = TT^(Xt{Xt + dC,,) + Re\,Pu,c) , (20) 

I y^us I ^ ^ 

with Xq = V*,Vgd, Xt = 1.464 ± 0.041, Pu,c = (0.2248)"^ Pu,c [33] and Pu,c = 0.41 ± 0.04 [34-36]. 



3 Observables 



The observables we are interested in, either to derive bounds from the existing measurements or to 
obtain predictions for future measurements, are the differential decay distributions of -B — > Xsi~^i~ 
and B -^ {K*,K)e^i- decays, and the integrated rates of Bs,d -^ ^^^~, B -^ X^-f, B -^ K^*'^iyi> 
and K — > tti^j?. In the following we analyse the theoretical expressions of these observables in 
terms of the (non-standard) Wilson coefficients of the MFV effective theory. 



3.1 B ^ Xsi+e- 

In order to minimise the theoretical uncertainty, we normalise all the observables to the cor- 
responding SM predictions. The numerical values for the SM predictions are computed at full 
NNLO accuracy, applying also QED corrections, following the analysis of Ref. [21,22], and non- 
perturbative 1/mi, corrections [19,23]. On the other hand, since full NNLO expressions for the 
non-standard operator basis in (7) are not available, the relative deviations from the SM are com- 
puted at NLO accuracy (with partial inclusion of NNLO corrections, as discussed below). Given 
the overall good agreement between data and SM predictions, and the corresponding smallness of 
non-standard effects, this procedure minimise the theoretical uncertainty. 

Following the standard notations of i? — > Xgi'^l" analyses within the SM, we define the 
effective coefficients 



Cf 



Cf 






a 



1 + ^cJ7(s) 



vr 
a,9 



Av + 



1 + —L09{s) [Ag + Tghimi, s) + Ugh{l, s) + Wgh{0, §)] + 

vr 



a.c 



IH UJg{s) 

vr 



^10 + • ^ 



(21a) 
(21b) 
(21c) 



where all quantities with a hat are normalized to the b quark pole mass (£ = x/ttij,). The 
leading-order gluon bremmsstrahlung corrections Wj(s) and the loop function h{x, y) can be found 
in [18,19], while the dots denote additional perturbative ag and electro-magnetic corrections, as 
well as non-perturbative power corrections, which we only consider for the SM normalization. 
The Wilson coefficients are evolved down to the low scale ^u^ = 2.5 GeV, using the input values 
in Table 2. As shown in [19], using such a low renormalisation scale minimise the NNLO QCD 
corrections to the B — > Xsi'^(.~ differential rate. An important difference with respect to the SM 
case is the presence of the Pq operator. To the level of accuracy we are working at, its inclusion 
is straightforward, since it is renormalised only in a multiplicative way. 

We are now ready to compute the inclusive decay spectra. Neglecting the strange quark mass 
while keeping the full dependence on the lepton mass, we get for the unnormalized forward- 



backward assymmetry and the rate 



i^FB 



ds 



ds 



256^5 \^tbVts\ {l-s} (1 



-2Re ( CfCf^* )s-ARe( CfCf^* ) + Re ( CfC^ ) rhg + 2i?e ( CfC^Q ' j m^ , 

(22) 



^^«-^V*.^;j^ (1-5)^/1-'^' 




+ (|C9^°|^ + |C 



«ff|2N 



+ 6m,^(|Cr|^-|Ci^o"r) 




+ Sm.ReiC'oC^^*) } ■ 



(23) 



As anticipated, these expressions are used only to evaluate the relative impact of new physics. Al- 
together, the estimated theoretical errors for the integrated inclusive rates reported in Table 1 are 
around 10% — 30%. In the numerical analysis this error is added in quadrature to the experimental 
one, which, at present, provides the largely dominant source of uncertainty.^ 



3.2 B^ {K*,K)i+r, B^ {K*,K)uu 

Exclusive modes are affected by a substantially larger theoretical uncertainty with respect to 
the inclusive ones. However, they are experimentally easier and allow to probe or constrain 
combinations of Wilson coefficients which are not accessibe (or hardly accessible) in the inclusive 
modes. The most interesting observables for our analysis are the forward-backward asymmetry 
(FBA) in S — > K*£~^£~ and the lepton universality ratio in i? ^> K^*'i~^£~ . In addtion, interesting 
predictions can be derived for the B -^ {K*,K)uD decay rates. 

In the B -^ K*l~^i~ case we need to introduce seven independent hadronic form factors {V, 



** There are two sources of theoretical uncertainties which, at present, are difflcult to estimate: i) non-perturbative 
power corrections of order asA/rrib', ii) radiative corrections associated to the soft-photon experimental cuts applied 
in the experiments. In the first case the effect is estimated to be of 0(5%) by simple dimensional counting [46] 
(more sophisticated estimates using the vacuum insertion method in the context of B — > Xs^ confirm this naive 
estimate [47]). As far as radiative corrections are concerned, the theoretical predictions include a correction due to 
coUinear logs [22] which does not correspond to the treatment of soft and coUinear photons performed so far by the 
experiments [48]. We guesstimate this mismatch as a 10% uncertainty in the high-g^ region and 5% error in the 
low-q'^ region. Given the large experimental uncertainties, which are still dominant, the impact of these additional 
theory errors turns out to be negligible. However, we stress that in the future a reduction of the theory uncertainties 
can be obtained only with a proper treatment of soft-photon corrections. 



Ao_2, ri__3 [31,49]). Following Refs. [50,51] we define the following set of auxiliary variables: 
2 



A{s) 
B{s) 
C{s) 

D{s) 

E{s) 
F{s) 
G{s) 

H{s) 

X{s) 



1 + rfiK* 

{l + rhK* 

1 



-ff^""^-^^^ , ^"^^.CfT^is) , 



Cf{s)V{s) + 



2m, 



■effr 



CT{s)Ai{s) + -,^(1 - rnK*)CrT2{s 



m 



K* 

■off. 



cff 



(1 - mK')Cf'{s)A2{s) + 2rh,Cf Ti{s) + 



m\ 



-Tois, 



leffr 



Cf{s) ((1 + mK*)Ai{s) - (1 - mK*)A2{s) - 2mK'Ao{s)) - 2mbCfn{s) 
2 



1 + mK* 

(l + mx*)Qt^i(5), 
1 



1 + mK* 



Ct^A^is) 



^C7it[(l + m^*)^i(s) 
mj, 



'1 -m/^*)^2(s) -2rfi/^*^o('S)] , 



(24) 



where now all quantities with a hat are normalized to the physical meson mass (and bremm- 
strahulng corrections are omitted from the effective Wilson coefficients without the tilde). Also 
in this case the above expressions are used only to evaluate deviations from the SM. As far as 
the SM predictions are concerned, in this case the most accurate results are obtained by means of 
QCD factorisation and SCET [28-30,52,53]. However, here we adopt a conservative point of view 
and evaluate the SM predictions and the corresponding errors by means of the parametrisation 
and QCD SR calculations of Ref. [31]. Then we enlarge the theoretical errors by the differences 
between our SM predictions and known QCD factorization and SCET results [28-30]. In practice, 
given the large experimental errors, this choice has almost no influence on the constraints derived. 
Using the auxiliary variables in Eq. (24), the un- normalised FBA then reads 



ds 



where 



Gp a^^ m^ |T/*T/ |2 ,-,/'s\2 ^ 

^yO^ \Vtsytb\ U{S) X 

sRe{BE*) + sRe{AF*) - 



-2 



m 



K- 



[XRe{XC*) - (1 - ?n^. - s)Re{XB*)] \ , (25) 



u{s) 




A = A(l, m^, , s) = 1 + m^, + s — 2s — 2m^, (1 + s) . 



(26) 

(27) 



Similarly the decay rate can be written as 

Ws y4' X { ^A {§ + mj) + ^A {s - ml) 






—2- W (A - u{sf/^ + 8m\4s + 2mj)) + jFp (A - u{sf/3 + 8m|^.(s - 4mj))] 



''K 

A 



"^K* 



[s|X|2 + 2miX{Re{XF*) - sRe{XH*) - (1 - mj^,)Re{XG*))] 



+ -^ [|Cp(A - iiisf/S) + |G|2 (A - u{sf/3 + 4m|(2 + 2m^. - s))] 
V [R^BC* + FG*)(1 - m|^. - s)(A - u(s)V3) + 4m^i?e(FG*)A] 



2rhj^ 



2 



2-^A [Re{FH*) - Re{GH*){l - m|.)] + jii'p-r^sA ^ . (28) 



In the B — > Ki^i~ case we only need three form factors (/+, /o and /t) [54]. Again we define a 
set of auxiliary variables [50, 51] 

A' = Cf{s)Ms) + -^CfMs), 

1 + mK 



D' = ^-^CfoiUs) - Ms)], 

X' = ^-^C'oW), (29) 



in terms of which 



A-^ 



ds~ 2i07r5 \^ibVts\ u{s)<{\A\ +|0|j 
+ 4|C'|2mf (2 + 2m\ - s) + 8Re{G'D'*)mj{l - m\) + A\D'\'^mis + |X'|2(2s - Am}) 
+ ARe{D' X'*)miS + ARe{G' X'*)rhi{l - m\) I . (30) 

In comparison to i? ^ K^*'£~^£~, theoretical calculations of the decay amplitudes for B -^ 
K^*'vv are considerably more reliable, owing to the absence of long-distance interactions that 
affect charged-lepton channels. We consider the missing energy distribution of the decay rates in 
terms the energy of the neutrino pair in the B rest frame and define the dimensionless variable 
X = Emissj which varies in the range (1 — m?,(,) )/2 < x < 1 — mj^(t) . For the B -^ K*^!? channel. 



10 



the partial rate then reads [55,56] 



dr 

dx 



\^t,bVtJ 



6C,p + Xt 



+s 



1287r5 

{rhK* + 1) Ai (s) + 






x'/' 



{rhK* + 1) Ai {§) 



niK* + 1 



niK* + 1 
{rhK* + 1)^ (^i^A + Am\, - 2rh\?j Ai (s) - A^2 {t 



2mK* {rhK* + 1) 



(31) 



while in the case of B 



Kuv we have 



d^- 128vr5 l^*'^^*'^l 



SC^p + Xt 



sm f^vi/ 



Ms)x'^' 



(32) 



3.3 S ^ Xsj 



Since the scalar-density operator does not contribute to i? — > Xgj, in this case the treatment is 
completely equivalent to the SM. We can thus take advantage of the complete NNLO analysis of 
Ref. [11]. Expressing the branching ratio in terms of the initial conditions of Cj and expanding 
around the SM value we can write [11, 12,57] 



B{B -^ Xsj)[E^yi.6 GeV] 



VtsV,l 



Vcb 



6.00 X 10" 



[3.15 - 7.186C7 + 4.7A{6C7) 



(33) 



where C^c = lKife/^fe|^r(-S — > Xc.ei')/T{B -^ XuSv) and its numerical value is reported in Table 3. 
Following Ref. [11] we assign a theoretical error of 0.23 x lO^'* to this expression. 

Thanks to the precise experimental measurement of the B — > ^^7 rate, Eq. (33) provide a very 
stringent bound on non-standard contributions to electric-dipole and chromomagnetic operators. 
As anticipated, we do not treat SCs as an independent parameter in the fit: the bounds on 6Ct 
derived by means of Eq. (33) should be interpreted as bounds on the linear dCj-dCg combination 
in Eq. (11). 

3.4 S,,rf ^ £+r 

The pure leptonic decays Bg — > i~^i~ receive contributions only from the effective operators 
Pio and Pq. These are free from the contamination of four-quark operators, which makes the 
generalization to the b ^ d case straightforward. 
The Bs — > i'^i~ rates can be written as 



where 



^^^^3 mBm^fBs 



m-. 



m 



Bs 



\1 + Ss\' 



+ 1 



m? 



m. 



\Ss 



Ss 



C^ml^ 



C^ m|^tan/?2 



mb 



AAiQ{mb + ms)mi XeXj, 4Aio{Huy rrib + rus 



(34) 



(35) 
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Quantity Value 



\Vus\ 


0.2255 ± 0.0007 [33] 


\Vcb\ 


(4.1 ±0.1) X 10-2 [10] 


\Vub\ 


(3.8 ± 0.4) X 10-3 [10] 


7 


(70.3 ± 6.3)° [9] 


IBs 


(0.260 ± 0.030) GeV [26] 


fBjfBd 


(1.21 ±0.06) [59] 



B{B^ XJu) (10.75 ±0.16) x 10^^ [10] 
Cue 0.58 ±0.01 [11] 

Table 3: Main inputs used in the numerical analysis (inputs not explicitly indicated in this Table 
are taken from Ref. [58]). 



The Bd — > £+£~ rates are obtained from Eq. (34) with the exchange {Vts, rns^, rng} — > {Vtd, itib^, 
rrid}- Neglecting tiny corrections of 0{ms/miy), this leads to one of the most stringent tests of the 
MFV scenario, both at small and large tan/3 values, namely the relation 



T(B, ^ i+n ^^ fB^niB,, 



r{Bd ^ £+i-) fB.mB, 



Vt 



ts 



Vt, 



(36) 



On the other hand, we stress that the relation between i?^,^ -^ ^^f- rates and IS.F = 2 observables 
discussed in [38] holds only in specific (constrained) versions of the MFV scenario. 

3.5 K -^ 71UU 

The branching ratio of charged and neutral K -^ ttvu decays can be simply expressed as 

B{K+ ^Ti+vu)=4{l + ^EM)\yA'' , B{KL^7T%u)=Kl[Im{y,)f , (37) 

where from [33,37] «;+ = 0.7867(43) x 10"^ and k^ = 3.3624(264) x 10"^ whereas the electromag- 
netic corrections A.em = —(0.30 ± 0.05)% . Similarly to Bg^ — > £"'"£" rates, also in this case we 
have two observables controlled by a single free parameter: the real coefficient dC^y in Eq. (20). 
Thus also in this case the ratio of the two K -^ ttvv rates leads to a significant model-independent 
test of the MFV hypothesis. 

4 Numerical analysis 

In order to determine the presently allowed range of the 5Ci we have performed a global fit of 
the AF = 1 observables in Table 1. The main numerical inputs beside rare processes used in the 
fits are reported in Table 3. The latter have been assumed to be not correlated and not polluted 
by new physics effects. In particular, as far the CKM angle 7 is concerned, we have used the 
results of Ref. [9] where the CKM matrix is determined using only tree-level observables. Our 
fitting procedure follows the method adopted by the UTFit collaboration [9]: we integrate over 
the probability distributions of inputs and conditional probability distributions of observables 
assuming validity of MFV to obtain (after proper normalization) the probability distributions for 
the 5Ci. 

The resulting ranges for the 5Ci and the corresponding bounds on the scale of new physics for 
the various operators are shown in Tables 4 and 5, and in Fig. 1. The bounds in Table 4 are the 
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5Ci 95% probability bound Observables 



6C7 


-0.14,0.06] U [1.42,1.62] 


B -^ Xsj, B -^ XJ+£- 


6C, 


-2.8,0.8] 


B -^ Xs£+i- 


6C10 


-0.4,2.3] 


B ^Xsl-^l", Bs^ fi+fi- 


5Cg/mb 


-0.09,0.09]/(4.2GeV) 


Bs -^ IJL^I1~ 


SCup 


-6.1,2.0] 


K^ -^ TT^UV 



Table 4: Combined bounds on the effective Wilson coefficients in the MFV scenario, and observ- 
ables used to derive the bounds (in the case of the scalar operator we report the bound in terms 
of the scale- independent combination SCg/nii,). 

Operator Aj@95% prob. [TeV] Observables 

Hi {DRXdXFC<Tf^uQL) {eF^u) 6l 

H^ {DrXAfc^^uT'^Ql) {9sGl,) 3.4 

(OlAfcTmQl) {eD^F,,,) 1.5 

i {QlXfc1,.Ql) h\jD^Hu 1.1« 



Xsl^B 


^ xj+i- 


Xsl,B 


^ Xsl+l- 


Xsi+r 




Xsi+r 


Bs -^ ^+/i 


xj+r 


Bs -^ ^+/i 


Xsi+£- 


Bs -^ ^+/i 


x,i+£- 


Bs -^ fJ.'^n 


XJ+£- 


B. ^ u+u 



i (OlAfct'^TmQl) hIt'^D^Hu 1.1 

{QL><FCl^lQL) {LLjf,LL) 1.7 

{QlXfcI^t'^Ql) {IlI^^t-Ll) 1.7 

(QlAfcTmQl) {En-i^ER) 2.7 

"A discrete ambiguity is removed at 90% probability, improving the bound to 2.3 TeV. 
Table 5: Individual bounds on the scale of new physics for the most relevant MFV operators. 

results of the global fit, where all the 5Ci are allowed to vary. The most interesting correlations 
among pairs of 5Ci of the global fit are shown in Fig. 1. On the other hand, the bounds in Table 5 
correspond to the bound on the scale of each non-standard operator, assuming the others to have 
a negligible impact. Note that the correlations of the 6Ci play a non-trivial role also in Table 5, 
by means of Eq. (6): each bound corresponds to setting one of the the a^ = ±1 and the others to 
zero. In case of sign ambiguities, the bound on the scale corresponds to the lower allowed value. 
In the case of the scalar-density operator, the translation of the bound on 5Cg into a bound 
on the scale is not straightforward as for the other operators. Assuming that the coefficient of 
051 in Eq. (1) does not depend on tan/? and setting 051 = ±1 we get 

T^^ = (1-5 TeV) X ( ^] [95% prob.] . (38) 

S max / ^ ^ 



max 



This bound, comparable to most of the bounds in Table 5, is especially interesting in specific 
models, where it can be identified with a bound on the mass of heavy Higgs fields. This happens 
for instance in the MSSM, where 051 is suppressed by l/167r^ (Csi being forbidden at the tree 
level by the Peccei-Quinn symmetry) but grows linearly with tan/3 [1]. In particular, setting 
|a5i|/A2 = tan/3/(167r2Af|f), leads to 

'tan/3^ ' 

"so" 



Mh > (830 GeV) x ( -^ J [95% prob.] . (39) 



As far the sign of C7 is concerned, we find that the wrong sign solution to C7 is still allowed, 
but has a lower probability compared to the SM sign. As shown in Fig. 2, the large contribution 
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ft 



I 68% prob. 
I 95% prob. 




prob. 0.15- 



B 95% prol 



0.5 



Figure 1: Correlation plots for the pairs of shifts in the Wilson coefficients. Allowed regions at 
68% (in green) and 95% (in red) probability. 




I 68% prob. 
■I 95% prob. 



i^ 



0.5 1 

6C7 



1.5 



Figure 2: Probability distribution of 6C7. 



to 6C7, corresponding to a sign flip of CY{mb), has a probability of about 30%. We stress that 
the sign flip of Cjlnif,) occurs only if Cg receives a sizable non-standard contribution. This is 
consistent with the conclusion of Ref . [60] , where the wrong sign solution to C7 has been excluded 
assuming small new-physics effects in the other Wilson coefficients. 

The impact of the low- and high-energy regions in B ^ Xsi~^i~ , which are often neglected, 
can be seen in Fig. 3, where we plot the most interesting 68% and 95% allowed regions with or 
without the information of these two measurements. In view of future experimental improvements, 
we report below the numerical values of the main observables expanded in powers of the 6Ci: 



B{B -^ Xs-i){E^ > 1.6GeV) = 3.13(23) x 10"^ x (1 - 2.286C7 + 1.516C^] 



(40) 
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region [1,6] GeV^ 
all regions 



SM 



3 




SCi 



-1 
5C9 



Figure 3: Correlation plots for the pairs of shifts in the Wilson coefficients, which are most affected 
by including the low and high energy regions in i? — > Xsi'^l~ ■ The left plot shows regions allowed 
at 68% (in green) while in the right plot, the regions shown are allowed at 95% (in red). The dark 
and light shaded regions correspond to the fit with and without the low and high energy regions 
m. B ^f Xsi~^i~ respectively. 



0.8(2) X 10~^ 
- 2.23SCr + O.235C9 



0.375^10 + 0M6C7SC9 



B{B -^ Xsl^r){q^ G [0.0021,0.04] GeV^) = 

X [1 + 2.025C| + 0.20(5C^ + Q.2{)5CIq 

+0.01(5C7^)2 - 0.01<5C7io5C^ + 0.01(5C7^] 
B{B -^ Xsi^£-){q^ G [1.0,6.0] GeV^) = 1.6(1) x lO"*' 

X [1 + 0.645C| + 0.615C| + O.GUCfo - 0.03607 + 0.866Cg 

+0.13{6C^f - 0.036Cio6C^ + 0M6C^] , 
13{B -^ Xs(.'^r){q^ G [14.4,25.0] GeV^) = 2.3(7) x 10"^ 

X [1 + 0.055C| + 0.285C7 + O.mSCl + 0.595(71^0 + 1-045C9 - 1.135Cio + 0.325C76C9 

+0.27{5C'i.f - 0.02(5Cio<5C^ + 0.02(5C7^] , (41) 



l.lbdCio + O.nSCjSCg 



B{Bs 



fi' fi 



4.1(8) X 10"^ X [(1 - 1.045Cio 



29.35C^f + 860{6C^f] 



(42) 



B{K+ -^ T^+vv) = 8.6(9) X 10"^^ [l + 0.96<5C^p + 0.2A5Cl^ 



(43) 



Since the experimental information of exclusive B — > K*l'^i.~ modes is quite good [27], we have 
also investigated the impact of including these observables in the fit. In this case we have used 
the results of Ref. [31] for the hadronic form factors and their corresponding errors, taking into 
account additional theoretical uncertainty due to neglected additional long distance effects [28-30] . 

The FB asymmetry plays a significant role. In particular, the normalised FB asymmetry 



AYB{q^ 



1 



d^FB 



dr/dg2 dg2 



(44) 
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Constrained fit 
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0.15 












^1 Constrained fit 
^1 Including ApB 




0.05 


^ 





^H ^B 


0.05 


^^ IH 


-0.1 




0.15 





-3 -2 -1 

(5Cn 



1 



Figure 4: Correlation plots for the pairs of shifts in the Wilson coefficients, which are most affected 
by including the exclusive Aps observables. Allowed regions at 68% with (dark green) and without 
(light green) exclusive Aps observables are shown. 



is interesting since the main uncertainty due to the form factors, namely the overall normalization 
of the decay rate, partly cancels out. Moreover, Afb{B -^ K*t^i~) in the low q^ energy region 
is very small in the SM due to the destructive interference of Cy and Cg (resulting in the famous 
zero of the asymmetry at low q^). It is therefore a good probe of the relative sign of the Wilson 
coefficients. The impact of including the presently available data on Ayb{B — > K*(.'^i~) is shown 
in Fig. 4.^ As can be seen, at present the additional information significantly reduces part of the 
ambiguities in the 5Cg-5CiQ and 5Ciq-6Cq planes only at 68% probability level. On the other 
hand, the overall bounds on the scales of individual coefficients do not change in appreciable 
way. This may seem at odds with conclusions reached in refs. [27,61]. However, the central 
experimental values for the high q^ region lie approximately 1.5 standard deviations above the 
range of possible theory predicitions within MFV satisfying other bounds. In addition, this range 
of theory predictions spans less than two experimental standard deviations. Under the assumption 
of MFV validity, the present FB asymmetry measurements therefore cannot significantly affect 
the 95% probability regions of 5Ci. As seen on Fig. 4, the situation would however improve 
dramatically, once the experimental precision would approximately double. 

Contrary to the FB asymmetry, it turns out that the K* longitudinal polarization is not very 
sensitive to new physics in the MFV scenario. 

For completeness, we report below the numerical expressions oi Afb{B — > K*i'^l~), integrated 



^Afb{B -^ K*£'^£~) has been measured both by Belle [61] and Babar [27]. However, in Ref. [61] only the fully 
integrated asymmetry has been reported. The normalised values of Afb{B -^ K*£^£^) in different q^ bins, as 
reported in Ref. [27] , represent the most useful information for our purpose. For this reason, we have restricted our 
numerical analysis only to the results in Ref. [27]. 
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Figure 5: 95% probability allowed ranges for the lepton universality ratio (muon/electron) for the 
differential distributions of B ^ Ki+i' (left) and B -^ KH+i' (right). 



and normalised to the decay rate as in Ref. [27] expanded in powers of the 6Ci: 



Afb{B ^ K*i+n,^^6.25 GeV^ 



-0.01(2) X 



T{B ^ A-^+^-)^2<6.25 GcV^ 

(1 - 20.^C7 - II..5C9 + l.^Cio + 21.6Cy6Cio + ll.6Cg5Cio) 

^ (1 + 1.26C^ - O.7SC7 + 0.5dCl + 0.55C72g + O.Q5C9 - l.OSCw + O.SdCjdCg) ' 



Afb{B -^ K*i+t 



Jq2>10.2A GeV^ 



0.20(8) X 



T{B ^ i^*£+£-),2>io.24 GeV^ 

(1 + O.SJCy + 1.25Cg - I.O6C10 - O.^dCjSCio - l.36Cg6Cio) 
^ (1 + 0.15C72 + OASCj + 0.66CI + OMC^q + 1.0<5C9 - 1.2(5C7io + OASCjdCg) 



(45) 



In the above expressions, we have neglected the scalar operator contributions, which are made 
negligible by the strong bound coming from Bg -^ fi'^ ^~ . However, 5Cg dominates possible NP 
effects in the lepton universality ratios of B ^ K^*'£~^i~ , as we discuss in the next section. 



5 Predictions 

Using the bounds on the MFV operators discussed in the previous section, we obtain a series of 
constraints for FCNC processes which are not well measured yet. The most interesting predictions 
can be summarised as follows: 

• R^^'/''\B ^ K^*k+£~). 

As pointed out in Ref. [51] the lepton universality ratios (or the muon to electron ratios of 
branching ratios) in i? — > K^*>l'^l~ are very clean probes of possible scalar density operators. 
Within the MFV framework there is a one-to-one correspondence between possible deviations 
from the SM in the lepton universality ratios and the contribution of the scalar-density 
operator in Bg -^ fj,'^fi~. Thanks to the substantial improvement on Bg -^ fi~^fx~, the effect 
is highly constrained. 

The present allowed ranges for the differential distributions of the lepton universality ratios 
are shown in Fig. 5. As can be seen, the effect is negligibly small in the B — > K* case. In 
the B ^ K case there is still some room for deviations from the SM, but only in the high-g^ 
region, where the decay is suppressed. The maximal integrated effect is reported in Table 6. 
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R'^P(Bs^yuV-) 



Figure 6: Correlations between Bg — > fJ,~^iJ-^, B^i -^ /^^^^, and Bg — > XgT^T". For a given 
value of B{Bs — > fi~^^~), normalised to the SM prediction (horizonatal axis), we plot the range 
of 13{B(i -^ ^+/i~) and 13{Bg — > Xs//+//~), both normalised to the correponding SM predictions 
(vertical axis). 



Note the major improvement with respect to the analysis of Ref. [51], where the loose bound 
on Bg -^ IJ'^l-i'~ allowed much larger non-standard effects. 

B{Bd ^ /i+//-) 

As anticipated, probably the best test of the MFV scenario in the AF = 1 sector is the ratio 

in Eq. (36). Taking into account the present bound on Bg -^ l-i-^IJ'^ this implies the upper 

limit 



B(B, 



fi^i-r) < 1.2 X 10-^ 



(46) 



at 95% probability. The clear correlation between B{Bii 
the MFV framework is illustrated in Fig. 6. 



/i'^^ ) and B{Bg -^ ^'^ ji 



m 



B{B ^ XgT+T~) 

The bound on the scalar-current operator allows us to derive a non-trivial bound also on 
B{B — > XgT^T~), which at large tan/? is sensitive to scalar-current operators. The para- 
metrical dependence on the modified Wilson coefficient of the integrated branching ratio 



IS 



B{B -^ XgT+T- 



)q2e[UA,25.0] GeV^ 



1.6(5) X 10" 



x(l + 0M5C^ + 0.385C7 + 0.81(5C| + OAUCIq + lA15Cg 
-0.785Cio + OASSCjdCg - 0M6Cw6Cl + 0.526CI + 0.326Cy 



(47) 



Assuming 5Cl 



(m. 



-/m^)5Cg, in agreement with Eq. (6), and taking into account the 



allowed ranges of the Wilson coefficients in Table 4, we obtain the 95% probability bound 
reported in Table 6. 

dApB/dq^B ^ Xg£+i-) 

Taking into account the available constraints on C7^9^io, the resulting allowed range for the 
inclusive FB asymmetry is shown in Fig. 7. As can be seen, in this case there is still a large 
room for non-standard effects: this observable is one of the few examples of quantities which 
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Observable 



Experiment 



MFV bound SM prediction 



rW^){B ^Ki^ 



1 0.17 ±0.28 [62,63]" [-0.004,0.14] 0(10"'') [64] 



rW^){B -^ K*£+i-) - 1 0.37+^:^;] ± 0.09 [63] [-0.002, 0.01] ' 



BiBa ^ fi+fi-) 



< 1.8 X 10- 



[24] 



< 1.2 X 10" 



W_ 
1.3(3) X 10 



^W 



B{B^XsT+T-) 



< 5 X 10" 



B{B -^ Kuv) 



[65] 



< 0.4 X 10" 



B{B -^ K*vu) 



[65] 



B{Ki 



TT^VV] 



[66] 



< 9.4 X 10"^ 

< 2.9 X lO-^'J 



1.6(5) X 10'"^ 
(0.5 ±0.1) X 10~^ 



(1.2 ±0.3) X 10^^ 
2.9(5) X 10-1^ 



"Here we quote naive averages of the values obtained by the experiments and with symmetrized errors. 

Table 6: Predicted observables. All bounds are 95% probability limits. In the first two lines, 
the SM predictions refer to the kineniatical region where lepton-mass phase-space effects can be 
safely neglected. The experimental results on uu modes are not explicilty indicated since only 
90% C.L. limits are available. 



could exhibit large deviations from the SM even in the pessimistic MFV framework. On 
the other hand, present data exclude most confgurations with flipped sgn(C9/Cio) at 68% 
probability. 

B{B -^ K^*^uD) and B{Kl -^ tt^uu) 

From the K~^ -^ vv bound on bCyy we can bound the rates of all the i^z^-type FCNC 
transitions^. The most interesting predictions are reported in Table 6 (in the B{B — > K^*>vv) 
case we take into account form-factor uncertainties along the lines discussed in Section 3.2). 
For completeness, we give the numerical expressions of the K^ — > 'k^vv and B -^ K^*'vv 
branching ratios: 



B{B -^ Kvu) 
B{B -^ K*uD) 



2.9(5) X 10"" [1 ± l275Cyp + OATSC^^] , 
0.5(1) X 10"^ [1 ± 1.26Cyp + OASC^ ^ 
1.2(3) X 10"^ [1 ± 1.26Cyp + 0A6Cl 



VV \ 1 



(48) 



6 Conclusions 

The MFV hypothesis provides an efficient tool to analyse flavour-violating effects in extensions of 
the SM. The effective theory built on this symmetry and and symmetry-breaking hypothesis leads 
to: 1) a natural suppression of flavour violating effects, in agreement with present observations, 
even for for new physics in the TeV range; 2) a series of experimentally testable predictions, which 
could help to identify the underlying mechanism of flavour symmetry breaking. 

In this paper we have presented a general analysis of the MFV effective theory in the AF = 1 
sector. From the current stringent bounds on possible deviations from the SM in various i?-physics 
observables we have derived the bounds on the scale of new-physics reported in Table 5. As can be 
seen, these bounds are perfectly compatible with new dynamics in the TeV range. We recall that 
the in models where new particles contribute to FCNC processes only at the loop level, the bounds 
on the particle masses are one order of magnitude weaker with respect to the the bounds reported 



^These predictions are vaUd only in the limit where we can neglect operators with the Y^YuYjYa flavour structure. 
This condition is always fuUfiled at small/moderate tan/3 and even at large tan/? holds in most excplicit MFV 



scenarios. 
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Figure 7: 95% probability allowed range for the inclusive normalised FB asymmetry in i? — > 
Xgi'^t' (left) and B — > K*t^t~ (right). The dashed curve denotes the central value of the SM 
prediction. The other three lines indicate the central values obtained flipping sgn(C7/Cio) and/or 
sgn(C9/Cio). 



in Table 5: A ~ 47rmNp. Thus in weakly-interacting theories respcting the MFV hypothesis and 
with no tree-level FCNC, we could expect new particles well within the reach of the LHC. 

Using the bounds on the effective operators, taking into account the correlations implied by 
the observables measured so far, we have derived a series of predictions for future high-statistics 
studies of flavour physics. This has allowed us to identify observables which could still exhibit 
large deviations for the SM even under the pessimistic hypothesis of MFV. The most interesting 
ones are the rare decays Bg^ — > fi~^fi~, rare B and K decays with a neutrino pair in the final 
state, and the FB asymmetry in B ^ K*i~^i~ decays. The allowed parameter space for the latter 
is shown in Fig. 7. 

Using the current bounds on MFV operators we have also identified a series of stringent tests 
of this symmetry principle. The most interesting negative predictions are summarised in Table 6. 
If these predictions were falsified by future experiments, we could unambiguously conclude that 
there exist new flavour symmetry-breaking structures beyond the SM Yukawa couplings. The 
effective theory allows us to obtain also some positive predictions, namely correlations among 
different observables which could still exhibit a deviation form the SM. The most interesting of 
such positive predictions is the correlation between Bg -^ l-i-^ fJ-^ and B^ -^ l-i-^ fJ-^ implied by 
Eq. (36) and illustrated in Fig. 6. A clear evidence of physics beyond the SM in both channels, 
respecting this correlation, would provide a strong support of the MFV hypothesis. 
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